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Abstract 

The results of a beam test performed at the KEK PS in June 2005 are presented. Drift 
properties of an argon-isobutane mixture have been accurately measured and compared with 
predictions at magnetic fields between and 1 Tesla. The r.m.s. point resolution of a padrow 
is compared with simulation and with an analytical calculation. The fundamental limitations 
due to detector geometry and gas properties are reviewed and the measured performances of 
the detector are found to be close to this limit. A numerical application to the case of a Linear 
Collider TPC is presented. 
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1 Introduction 



Three of the four detector concepts which have been proposed for the Linear Colhder foresee a large 
Time Projection Chamber (TPC) as a main tracker. This allows continuous tracking to be performed, 
yet with a minimal amount of matter. Depending on the detailed designs, the TPC should have 
about 200 padrows with a space resolution between 100 and 150 microns in the R0 direction. Mainly 
three technologies are currently considered for the gas amplification in these TPCs : a Multi-Wire 
Proportionnal Chamber (MWPC), a Micromegas chamber, and a multiCEM structure. To this 
end, R&D has been pursued since the beginning of the decade within the LC-TPC collaboration [1] 
and has lead to the construction of several prototypes. It was felt useful to gather around a single 
experiment and to have a way of comparing the various technologies in a well-defined framework, 
with the same readout electronics and chamber geometry. This is the purpose of the Asian- Canadian- 
European Multi-Prototype collaboration. It started with the construction of a MWPC chamber in 
MPI Munich, which took beam data at KEK in April 2004. Then the chamber was equipped with a 
triple-GEM structure and took beam data in April 2005. This paper reports on a beam test carried 
out with a Micromegas endplate in June 2005. 

The detector and the operation conditions are described in Section 2, where an assessment of data 
quality is given. The results on gas properties (drift velocity, diffusion) are presented in Section 3. 
Section 4 is devoted to the spatial resolution of the device. It starts with a theory elaborated in our 
group which is then compared to simulations and measurements. Fundamental limitations from the 
gas mixture and pad geometry are assessed and consequences for a Linear Collider TPC are drawn. 

2 Experimental Setup and Data Taking 

This experiment was carried out using the 112 beam line at the KEK 12 GeV PS. The beam line 
provided a secondary beam of pions or protons with momenta up to 4 GeV/c through the interaction 
of 12 GeV protons on a Be target, followed by the charge and momentum selection with a set of 
dipole magnets. The beam spill has a flat top of 1.5 s with a repetition rate of 0.25 Hz. There were 
4 scintillation counters, TCI through TC4, whose 4-fold coincidence triggered the data acquisition 
on 4 GeV negatively charged pions. A typical trigger rate was 12 Hz. The first two trigger counters, 
TCI and TC2, were placed at the entrance of the beam just downstream of a beam slit to control 
the beam intensity and had an overlap region of 2 x 2 cm^. The other two trigger counters, TC3 and 
TC4, were located at 8 m downstream of TCI and TC2, that is, just in front of our Multi-Prototype 
TPC (MP-TPC). TC3 and TC4 had an overlap region of 30 x 10 cm^ that matched the drift region 
of the chamber. The MP-TPC together with TC3 and TC4^ were installed in a Persistent Current 
solenoidal Magnet (PCMAG) having a bore diameter of 85 cm and a length of 1.3 m with a very thin 
wall of 20 % radiation lengths. The magnet capable of creating a field up to 1.2 Tesla was operated 
in the closed loop mode and provided a field uniformity better than 0.5 % in the drift region of the 
MP-TPC that was ahgned with the magnet axis, so that the electric drift field was parallel to the 
magnetic field. 

In the following, unless otherwise stated, the beam was shot perpendicular to the drift axis of the 
MP-TPC. The nominal size of the beam in these conditions is about 4 cm at the chamber. To cover 
the whole active volume of the MP-TPC, we hence spread the beam by inserting a 5 cm-thick lead 
brick just downstream of TCI and TC2 in the normal data taking conditions. 

^TC3 and TC4 were equipped with fine-mesh photo multipHers, Hamamatsu R6682, which allowed operation in a 
magnetic field up to 1 Tesla. 
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Figure 1: The Multi-Prototype TPC (MP-TPC) 



2.1 Multi-prototype TPC 

The MP-TPC has a cyhndrical drift region of 261 mm in length and 145 mm in diameter. As seen 
in Fig. 1, the field cage is made of fourteen 15 mm-wide circular rings spaced 3.2 mm apart and a 
cathode disk with a small hole at its center^. All of these are embedded into a gas-tight vessel, which 
is closed on one side by a detachable endplate detector that carries a Micromegas foil and readout 
pads. The ring closest to the Micromegas detector is at a distance of 6 mm from the mesh. 

The Micromegas endplate, built in Saclay and Orsay, consists of a 0.8 mm thick Printed Circuit 
Board bearing anode pads, with a mesh stretched on a frame and mounted above the pad plane. The 
mesh is 5 /xm thick Cu having circular holes of 35 /im in diameter placed with a 60 /xm pitch. The 
100 X 100 mm^ frame leaves a 75 x 75 mm^ active area. A 50-micron gap is maintained between the 
mesh and the pad plane by kapton pillars. In normal operation a voltage of about 350 V is applied 
to the mesh relative to the pad plane at V. The resulting electrostatic force sticks the mesh to the 
anode plane. There are 12 rows of 32 pads on the anode plane. Each pad has a rectangular shape 
and is placed at a pitch of 6.3 mm along the beam [y direction) and 2.3 mm transverse to the beam 
[x direction). 

The chamber is filled with an argon mixture containing 5% isobutane and is operated at room 
temperature and at atmospheric pressure. 

2.2 Readout Electronics and DAQ 

The 12 rows of 32 readout pads (384 in total) are connected via 30 cm long flat cables to 24 ALEPH 
preamplifiers [2], each having 16 channels and reading out a half of a pad row. The shaped signals 

^The hole lets X-rays from an ^^Fe source enter the drift volume, for gain calibration and monitoring purposes. 
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Figure 2: A typical 4 GeV pion track from the KEK ^2 beam seen by the Micromegas MP-TPC in 
a 0.5 T magnetic field ; the squares represent the pads hit with a colour code corresponding to the 
charge deposited. 

from the preamplifiers are sent to ALEPH TPDs (TPC Digitizers) in a Fastbus crate via 15 m long 
twisted pair cables and are sampled at a rate of 12.5 MHz and digitized over 8 bits. The data are then 
read out, via a Fastbus- VSB translator unit, FVSBI9210, by a VME board computer, F/C8234, 
operating on OS/2. The readout data are stored via TCP/IP connection on a Linux PC in the LCIO 
format [3]. The data acquisition time is about 4 s per event, which limited the data acquisition speed 
to about 1 event per spill. 

2.3 Data Taking 

A tj^ical event taken at a magnetic field of 0.5 T is shown in Fig. 2. 

During the data taking, the gain was continuously monitored by a ^^Fe source. The mesh signal 
was readout by a fast charge amplifier ORTEC 142B and the signal was sampled with an AMPTEK 
Multi-Channel Analyser MCA8000. A spectrum obtained this way is shown in Fig. 3 where the 
5.9 keV line and the escape line in argon are seen. The source was not coUimated. The Landau 
distribution of the ionisation produced by the 4 GeV pions from the beam is also visible and peaks 
at 12.5 keV as expected for the 75 mm active length of gas. Over a period of 60 hours the extreme 
variations of gain were ±3% and the r.m.s. gain variations were 3 per mil over this period. The 
gain is measured as a function of the mesh voltage (Fig. 4). The data presented here were taken at 
a mesh voltage of 320 V, corresponding to a gain of 3650, in a magnetic field of 0.5 and 1 T, and a 
mesh voltage of 340 V (gain of 7500) for the data taken with no magnetic field. 

A total of 4020 triggers were collected at 5 = OT, 6111 at 0.5 T and 5166 at 5 = 1 T. The 
temperature was 28 to 32 °C during the data taking. The gas pressure in the vessel followed the 
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Figure 3: Energy distribution from the mesh signals from the source and from the beam 
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Figure 4: Gas gain as a function of the mesh voltage 
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atmospheric pressure which was stable during the runs used in this paper. 

In the following analyses, the first two and last two padrows. as well as the four leftmost and 
two rightmost pads of each row, are not used in the measurements, to avoid edge effects. However, 
when a hit is found close to the edges of this fiducial region, the neighbouring pad(s) in the region 
are recovered and used in the track reconstruction. 

The azimuthal angle distribution is shown in Fig. 5a. As stated earlier a wide beam was obtained 
by adding a 5 cm-thick lead brick in the beam line. The resulting z distribution of the tracks is 
shown in Fig. 5b. The uniformity of the detector can be assessed from the average residual vs track 
position in x (Fig. 5c) and vs padrow number Fig. 5d for the 8 fiducial padrows. Distortions up to 
50 microns are observed at non-zero magnetic field (likely to be due to E x B effects), but they do 
not affect the resolution measurements where only the r.m.s. of the residuals is considered. 

2.4 Data Analyses 

The exactly same data sample has been analysed with two independent analysis methods and com- 
pared in the following. They differ primarily by the way the track parameters are fitted to the pad 
data in the x-y plane. The geometrical track parameters are the intercept of the track with the y — 
axis (xq), the azimuth at this point (0o)) and the inverse radius of curvature (l/R). 

In the first method, a global maximum-likelihood fit of the geometrical track parameters, together 
with the track width (cxtrack) to the charge of the pads is carried out [4]. In other words, a track is 
represented by an arc of a circle with a gaussian charge spread of width Utrack along the x direction, 
and (Ttrack IS an additional fit parameter. The pad charge and time are determined as follows: the 
ADC counts of the 3 preceding and 3 following the time bucket where the maximum is reached are 
added together with the bucket of maximum charge. This sum is required to exceed a minimum 
of 7 ADC counts for the pad to be included in a hit. The time associated to the hit pad is the 
charge- weighted average of the seven buckets used in the charge integration. 

The log-likelihood function which is maximised with respect to the 4 track parameters reads: 
Snjln(/j), where is the fraction of the charge expected in pad i, obtained by integration over the 
pad of a function with a gaussian profile transverse to the track, rii is the number of electrons before 
amplification corresponding to the charge readout on the pad i. This method is implemented in the 
FORTRAN 95 program FTPC[5]. 

In the second method, a fit is performed to hit points on a row-by-row basis. The method 
is implemented in C++ and available as the Double-Fit program[6]. The Double-fit program starts 
with cluster finding in the plane of pad-number vs time-bucket for each pad row. To build a 2- 
dimensional {x and z) cluster, consecutive time buckets are summed, starting from the first pulse 
encountered above 9 ADC counts, up to the first below 6 counts. Such a hit is accepted if its peak 
pulse height exceeds 15 ADC counts. The coordinates of the hit is then calculated as its charge 
barycenter. Following the cluster finding, the program then performs a fit of either a straight line 
at -B = or a circle at B = 0.5 or 1 T. 

3 Gas Properties 

3.1 Drift Velocity Measurements 

Measuring the drift velocity of gas mixtures is an interesting test of the electron transport simulation 
in gases, or alternatively can be considered as a check of gas purity and composition. The drift field 
was set to 220 V/cm and the detector was taken outside the magnet, itself not energized. For this 
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measurement a novel and very simple technique was used: the beam was shot at an angle of 45 degrees 
to the center of the cathode. A 1 x 1 cm^ scintillator was added in the trigger to select pions crossing 
the cathode plane in the center. Electrons from ionisation close to the cathode drift all the way along 
the axis of the TPC reaching the central pads of the detector, in the region where no electric-field 
distortions are expected. The TPD measures the time elapsed between the trigger arrival and the 
signal arrival. The end-point of this time distribution, added to the delay between the trigger and 
the readout, is the time taken by the electrons to drift along the 260.8 mm chamber. The drift time 
is found to be 5, 907 ± 30 ns. The trigger delay was measured to be 310 ± 5 ns. Adding a rough 
estimate of 30 ± 20 ns for the NIM to ECL converter at the entrance of the TPD, and subtracting an 
estimated 75 ns for the pad signal delay in the 15 m readout cables, the total drift time is measured to 
be 6.172 ±0.045 /xs. Dividing the chamber length by this time yields I'drift = 4.226 ±0.031 cm///s. An 
alternative method has also been used, with a wide beam perpendicular to the axis of the chamber 
and also measuring the endpoint of the time distribution. It leads to the value of 4.157±0.036 cm/ /is, 
in agreement with the former. These two values are combined assuming fully correlated systematics 
to yield the measurement of the drift velocity of electrons in Ar+5% isobutane mixture at an electric 
field of 220 V/cm of 

Vdrift = 4.181 ± 0.031cm//xs. 

This value is used for the determination of the z coordinate of the hits along the axis of the chamber 
in the following data analyses. 

This measurement is in very good agreement with the Magboltz[7] prediction of 4.173±0.016 cm///,s. 
The uncertainty on the prediction is dominated by the error on the gas composition, stemming from 
a 2% uncertainty on the isobutane gas flow. No consideration on the model, the approximations nor 
on the input data used in the simulation enter this estimate of the prediction uncertainty. 

3.2 Transverse Diffusion Constant Measurements 

For large enough drift distances, so that the ionisation charge is spread over several pads, the width 
of the track crtrack allows the determination of the diffusion constant Cd using the relation (Jtrack = 

In the global likelihood method, Utrack is obtained by maximising the track likelihood simultane- 
ously for the geometrical track parameters. Fig. 6 shows for the 3 values of the magnetic field the 
relation between the average crtrack the drift distance z. The expected linear dependence on the 
drift distance z. for large z^ is clearly seen. The data at low z. where the sensitivity to Ctrack is lost 
and where the spread in track widths is comparable to the average track width, are not included into 
the hnear fit. The shght offset in the case of zero magnetic field, corresponding to 0.3 mm added in 
quadrature to the track width, can probably be attributed to delta rays. At higher fields the path of 
these delta-ray electrons are expected to be limited by the magnetic field. The fitted Cd values are 
given in the first line of Table 1. The quoted uncertainties include systematics from the noise level 
and threshold used to reconstruct the hits. 

In the second method, the average fraction of the charge falling on a given pad is plotted as a 
function of the distance between the track and the pad center. It can be shown (see Appendix A) that 
for large enough drift distances and for high enough magnetic fields to curl up delta-ray electrons, 
this distribution becomes gaussian and its width scales as gpr = ^Jw'^ /12 + C^, ^, where w is the pad 
pitch. The measured values of C/j by this method are shown in the second line of the table. These 
two measurements are in good agreement with each other and with the Magboltz prediction, at the 
level of a few percent. 
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drift distance (mm) 

Figure 6: Track width as a function of the drift distance for B = 0, 0.5, and 1 T 

Table 1: Diffusion constants in fim/y/cm measured by two methods and Magboltz prediction 



Magnetic field 


T 


0.5 T 


1 T 


Global likelihood 


488 ± 11 


314 ± 15 


209 ±7 


method 


475 ±5 


293 ±5 


194 ± 18 


Magboltz 


469.3 


284.1 


192.6 



4 Resolution 

4.1 Theory and Simulation 

4.1.1 Effective Number of lonisation Electrons 

In this section, we present an analytic formula for the spatial resolution of a TPC equipped with 
a readout plane consisting of a gaseous detector that amplifies track electrons and rows of readout 
pads to measure their charge centroids. For simplicity, consider a charged particle passing through 
the TPC at the right angle to a pad row and at a drift distance z from the readout plane. In this zero 
crossing angle case, as far as the size of the primary ionisation clusters is negligible^, all of the track 
electrons created in the ionisation will have a 5-function-like distribution peaking at, say, x = x, 
when projected onto the x-axis, namely, the axis in the pad row direction. 

These track electrons drift toward the pad row. While drifting, each of these track electrons 
experiences transverse diffusion and will have an a;-coordinate deviated from x = x hj Axi, i.e. 

•^This assumption is justifiable at high magnetic field which is expected to curl up delta rays and hence suppress 
the size of the primary ionisation clusters. 
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Xi — X -\- Axj, according to the probability distribution: 

where the subscript, i, means i-th electron and Ud — Cd^/z with Cd being the transverse diffusion 
constant. 

The number of track electrons that will reach the gas amplification region in front of the pad row 
fluctuates statistically. Let us denote the probabihty of getting N such ionisation electrons eventually 
contributing to the signal induced on the pad row by Pi{N; N) with N being the average: N — (N). 

Each of these A'^ ionisation electrons will be gas-amplified at the readout plane by a factor, G, 
which is assumed here to fluctuate according to a Polya distribution: 

with G being the average gas gain: G = (G) . Notice that the Polya distribution becomes exponential 
in the ^ — > limit, while it coincides a 5-f unction in the ^ — > oo limit. Notice also that the Polya 
distribution has a variance: cr^yg = 1/(1 + 9), which goes to zero as ^ — > oo as expected. 

For illustration purpose, let us assume, for a while, an idealistic readout plane which measures 
the x-locations of individual electrons with infinite accuracy but with relative weights of gas gain 
values. Then, the center of gravity of these N electrons at the readout plane will be given by 

~ ^ — X-{- . {O) 

The charge centroid, x, will then be distributed according to 

oo N , I. . 

P{x- x) ^ Pii^-^ ^) n / dAxiPniAxf, a^) / d{G,/G) PaiGjG; 

N=l 1=1 ^"^ 

Eili G,A x. 
Ef=i G 



><M^-^-^^^ir^l- (4) 



Under the assumption that is large enough and hence 

^-^E^. (5) 



1=1 



we obtain the variance of the center of gravity, cr^, by inserting Eqs.(3) and (4) into its definition 
and carrying out the integral: 
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J dx P{x; x) {x — xY 



where use has been made of 



^■^^ I N=l 
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and 

((§)} = / d{G/G)PG{G/G-e) (8) 
The A^eflf is hence given by 




Notice that A^efr is in general significantly smaller than TV = {N) due to ionisation statistics. 
The gain fluctuation further reduces A^eff[9] by as much as a factor of two for the exponential gain 
fluctuation: 6* = 0. 



In order to find a lower limit on d, let us consider the total charge distribution. The total charge 
after gas amplification is given by 

N 

Q = Y.Gi (10) 
1=1 

and then its probability distribution function by 

Pq{Q) = E Pi{N-. N)I[(f d{G,/G) Pg{G,IG- d)\ Aq-Y^g\. (11) 
Prom this we can readily obtain the variation of the total charge as 

o% = JdQPQ(Q)Q^-[JdQPQ(Q)Q 

This imphes that in the case of a Landau-like Pi{N; TV), for which N <S the gas gain fiuctuation 
will not affect the total charge distribution very much. 
For X-rays from ^^Fe, however, we expect 
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^ ' ^ +F], (13) 



N \l + 9 

where the Fano factor F is about 0.2 and N is about 220 for argon. The best resolution so far attained 
with a Micromegas detector with exactly the same type of mesh is about 6.8% in r.m.s.[8], which 
implies 6 ^ 0.22. This is a lower limit in the sense that it neglects experimental effects increasing 
the width, other than gas gain ffuctuations (electronic noise for instance). However, 6 can be smaller 
if the Fano factor is less than 0.2. Taking a value F = 0.16, still consistent with measurements and 
theoretical calculations [10], would lead to a minimum value of 0.166 for 9. 
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4.1.2 Effects of Finite Size Pads 

In the last section, we considered gain fluctuation in the gas amplification process at the end plane 
detector but assumed that we could measure the location of a single electron with infinite accuracy. 
We now introduce a pad row of pitch w to measure the charge centroid: 

^ = i:qj{wj)/j:qv (14) 

where Qj is the charge on pad j and is given as the sum of contributions from N seed electrons: 

N 



i=l 



with fj being the response function of pad j for seed electron i arriving at the location x + Axi and 
AQj being the electronic noise on pad j. Notice that the pad response function is normalized as 

Y,f,{x + Ax,) = 1. (16) 
j 

The probability distribution for Qj is then given by 

PjiQr,^) = Y,Pi{N-N)Jl(f dAx,PD{Ax,-aa) f d{G,/G)PG{G,/G-e)] 

N=l i=l J / 

X J dAQjPE{AQj;aE)s(^Qj-Y.Gi-fj{x + Axi)-AQj^, (17) 



where Pe represents a constant electronic noise with (AQj) = and (AQj) = aj;. On the other 
hand, the probability distribution for the charge centroid is given by 



P{x- x) = Pii^-^ ^) n / dAxiPoiAxi-, aa) / d{G,/G) PoiGi/G; 9) 

Y[(^J dAQj Pe{AQj; aE) J dQj S (^Q, - E " fji^ + ^^i) ' ^Qj^ ) 



X 

j 



which replaces Eq.(4). 

Since the probability distribution P{x; x) depends on the true location of the seed cluster x, we 
average over x to define cr^: 



4 = 1^'^' d J dx P{x; x) {x - xf. (19) 



Substituting Eq.(18) in this, and performing integration over AQj, Axi, and AGi/G in this order, 
and then averaging over A^, we obtain 



X 



oi = r'^d[-\ X 



Y^^jw) {fj{x + Ax)) -X 
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where we have ignored the electronic noise as compared to the total charge: 

N N 
j i=l j j i=l 

and 

EG, ~ NG 

i=l 

as usual. 

Notice that the pad response function only appears in the following two forms: 

{fj{x + Ax)) = I dAxPoiAx; a^) fj{x + Ax) (21) 

and 

{fj{x + Ax)h{x + Ax)) = I dAxPoiAx; a^) fj{x + Ax) U{x + Ax), (22) 

and can be numerically evaluated, once the functional form of fj is given. 

The formula, Eq.(20), can be qualitatively interpreted as follows. The first term is the mean 
square of the difference between the charge centroid and the true location of the seed cluster and 
is independent of the number of primary electrons^. This term vanishes in the narrow pad limit, 
w — s> 0, while it approaches the famous {w / \/V2)'^ in the wide pad limit, w 3> (7^. The second 
can be interpreted as the combined effect of the diffusion and the gas gain ffuctuation. As we have 
seen in the previous sections, their contributions scale as {1 / N) (^{G / G)"^^ = 1/A^eff- The last term 
represents the contribution from the electronic noise and is independent of the shape of the pad 
response function or the diffusion. It scales as {w a e / G)"^ {1 / N'^) . 

4.1.3 Application to a Micromegas-like Readout Plane 

The qualitative observations we made in the last section agree with naive expectations. For quan- 
titative comparison with data, however, we need a concrete form of the pad response function. For 
simplicity, let us assume that the spatial size of the avalanche caused by a single seed electron is 
negligible compared to the pad width, as expected for a Micromegas-like readout plane. In this limit, 
the pad response function becomes hodoscope-like: 

fj{x + Ax) = e {{x + Ax)/w - J + 1/2) e (j + 1/2 - (x + Ax)/w) , (23) 

for which we have 

{f,{x + Ax)Ux + Ax)) = {f,{x + Ax)^)6,k= {f){x + Ax)) 6,^. (24) 

''It is well known that there is an S-shape systematics in the difference between the simple charge centroid and the 
true cluster location. This term can hence be eliminated by correcting the charge centroid for the 5-shape systematics. 
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Substituting these in Eq.(20), we obtain 

r+l/2 



^' - /_" 



dl-] X 
1/2 \W, 



JlU'^) {fj{x + Ax)) -X 



+ 



N 



I 



iw + Ax)) - ^{jw) {f,{x + Ax)) 



V 



(25) 



As long as ad ^ w as at long distance, in the integrand of the above equation the first term can be 
ignored and the second term can be approximated (see Appendix B) by 



Y^ijwf {f,{x + Ax)) - x' 



C7i + 



12 



(26) 



In this long drift distance limit, we have an asymptotic formula: 



4 - 



Q. 2N 



w 



12 



(27) 



which implies 



(7o = 



w 



w 



12, 



12/ VN^' 



if the electronic noise is negligible. 

The integration in Eq.(25) can be carried out numerically at any drift distance, using 



{fj{x + Ax)) = J dAxPD{Ax;ad)fj{x + Ax) 



rwj—x-\-w/2 ][ 

/ dAx 



exp 



' Ax'V 
2 



/(i + l/2V-x\ _ / (j - 1/2)^ - ~x 



(28) 



(29) 



This and Eq.(25) imply that o^jw is a function only of Gdjw and A/'eff- Fig- 7 plots a^jvo as a 
function of this scaling variable, ad/w. Notice that the full theory curve merges into the asymptotic 
formula at around Od/w ~ 0.4, which means that the effect of finite pad pitch becomes neghgible for 
(^d/w > 0.4. The full theory has a fixed point, a^/w — l/-\/l2, at Od/w = 0, while the asymptotic 
formula scales as l/^/N^s- The full theory curve attains its minimum of about cr^/w — 0.1 at around 
Gd/w ~ 0.3. 



4.1.4 Comparison with Monte Carlo Simulation 

A Monte-Carlo simulation has been carried out to check the analytical theory and to estimate the 
expected effective number of electrons for the gas mixture used in the tests. 
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Figure 7: Expected spatial resolution normalized by the pad pitch for A^eff = 18.5 as a function of 
the scaling variable ad/w = Coy/z/w. 

In a first step track segments are generated with a uniform distribution of x across a 2.3 mm 
pad and for a set of z values along the drift axis. The number of ionisation clusters is generated 
along a 6.3 mm pad according to a Poisson distribution with mean value 32.0. This average number 
of clusters is estimated as follows: a concentration-weighted average of the number of clusters from 
minimum ionizing particles (m.i.p.) in Ar (23 e~/cm) and in isobutane (84 e~/cm), taken from 
ref. [11], multiplied by the calculated ratio of ionisation dN/dx for 4 GeV pions with respect to 
m.i.p., taken to be 1.23 from ref. [12]. In each cluster a number of electrons is generated according 
to the argon cluster size distribution given in ref. [13]. The average 1/A^ is predicted to be 1/38.7. 

Then each individual electron is transported over the distance z. The transverse diffusion in the x 
direction is simulated by varying a; by a random amount, following a gaussian law, the width of which 
is given by C^y^, where Cr, is the Magboltz prediction for the diffusion constant : Co = 0.469, 0.285, 
and 0.193 mm/ ^^cm for 5 = 0, 0.5, and 1.0 T, respectively. 

The last step is to simulate the gas amplification gain. For this a Polya distribution with the 
parameter equal to 0.22 is used. For every track, the charge sum on each pad is calculated. Then 
hits are reconstructed as barycentres of the pads hit. The x resolution for each z is plotted in Fig. 8. 
The analytical formula, with input A^cfr = 21.3, as given by Eq. 9, reproduces the Monte Carlo data 




Asymptotic Formula 
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Figure 8: Expected spatial resolution with readout pads at a pitch of -u; = 2.3 mm for (l/N) = 1/38.7 
and 9 = 0.22, assuming Magboltz results Cd = 0.469,0.285, and 0.193 mm/y^cm for B = 0,0.5, and 
1.0 T, respectively. 



very well. 



4.2 Comparison with Measurements 

In the global likelihood analysis, the hit position on each row is estimated by a fit to the charge 
depositions in the row, with the track parameters fixed to their value using the 8 fiducial rows. The 
resolution is then calculated as the geometric mean of the r.m.s. values of the distributions of the 
residuals with and without the pad row in question in the fit (see Appendix C) so as to eliminate 
contributions from the tracking errors. To avoid being sensitive to outliers (produced by noise for 
instance) the hits situated at a distance of more than 4 standard deviations from the track are ignored. 
Alternatively, a gaussian function was fitted to the whole residual distribution to estimate its r.m.s.. 
This second estimate of the resolution is lower than the former by only 7%. The resolutions at 
-8 = 0, 0.5, and 1 T are shown as a function of the drift distance in Figs. 9a), b), and c), respectively. 
The results of the global likelihood method are the triangles. The method results are also shown 
(square data points). The two methods agree each other very well at long drift distances. The 
discrepancies at short distances will be discussed later. The likehhood method makes a better use of 
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Table 2: Effective number of electrons A^eff measured by the two methods 



Magnetic field 


T 


0.5 T 


1 T 


Global hkelihood 


16.5 ±3.4 


18.1 ± 1.2 


22.8 ±3.0 


fit + barycentre 


15.1 ± 1.2 


18.7 ±2.6 


15.7 ±7.7 



the pad signal information. 

These measured resolutions are compared with the theoretical predictions explained in the pre- 
vious section, with A^efr = 18.5. The fit using the barycentre of each row leads to values of A^eff 
consistent with this (see Table 2). The theory reproduces the data very well when Co^/z/w > 0.4 as 
expected for B = 0.5 and 1 T. The theory seems to underestimate the resolution even at longer dis- 
tances at B = and the discrepancy seems independent of the drift distance. This discrepancy could 
be attributed to the finite geometrical size of primary ionisation clusters that would be significant 
only when the curhng up effect of delta-ray electrons is absent. 

The hollow data points at small drift distance are not used in the fit. The reason to discard 
them is that they are biased toward low values, as the hits arc reconstructed towards the middle of 
the pads. An example of this is given in Fig. 10, where the distribution of the position of hits with 
respect to the middle of a pad is shown; it is clear from this plot that for too small drift distance the 
hits are preferentially reconstructed in the middle of a pad, so that the whole track is biased towards 
the middle of the pad and the residuals are underestimated. This is the main reason why the rise of 
the resolution at low z is underestimated in the data (this effect is less prominent in the case of 1 T 
where the track curvature in the magnetic field plays the role of an effective pad staggering). This 
effect is present in both methods and could be avoided only if we had had an external measurement 
of the track position. Notice, however, that unlike the method which uses the charge barycentres, 
the hkelihood method is free from the -S-shape systematics (the first term in Eq. 20) as long as each 
pad row has multiple pads above threshold. This is why the likelihood method tends to deviate from 
the theory and give better resolutions than the method in the short drift distance region. 

The fitted values of A'gff are given in Table 2. Considering that A^eff is independent of the magnetic 
field, the 3 data sets can be combined to obtain a measurement of A^eff — 18-5 ± 1.1. This value can 
be compared with the expectation A^eff = 21.3 ± 2.7 obtained in the previous section fixing the gain 
fiuctuation parameter 9 to 0.22, where the error is evaluated by varying all the other input parameters 
within their admissible range and adding in quadrature the variations of A^eff- The agreement is good, 
and favors significant gas gain fiuctuations. 

As explained in the previous section, a non-zero value of cto is expected from the fit to the data 
points of the functional form 

a = + Clz/N^s- 

The combined measurement of ctq is 165±18 /xm, in good agreement with the expectation w / \/VlN^ = 
154 iim. 

4.3 Extrapolation to ILC-TPC 

Conforted that the theory of resolution presented here is a pertinent approximation for a Micromegas 
TPC, we can use it to predict what can be expected in the 4 Tesla case, as proposed for the Linear 
CoUider TPC. Fig. 11 represents the expected resolution as a function of z, for two values of the 
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Figure 9: Resolutions at (a) 5 = 0, (b) 5 = 0.5, and {c) B = 1 T, as a function of the drift length 
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Figure 10: Distribution of the distance between the track and the center of a pad for a small drift 
distance sample (dots) and a large drift distance sample (triangles) from the global hkelihood analysis. 
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Figure 11: Resolution in the Linear Collider case as a function of the drift length, for two values of 
the pitch, for the gas used in this study 



pitch for the gas studied in this paper. One can see that at such a high magnetic field, the diffusion, 
with a constant of 63 jimj y^cm, is not sufficient to spread the charge enough. This results in the fact 
that the hodoscope effect is felt up to distances over one meter in the case of a 2.3 mm pitch. Even 
a 1 mm pitch, which would require a very dense readout electronics and would feature pads with a 
large aspect ratio, would not allow the target average resolution of 100 jJLva. to be reached. 

In the case of the triple mixture Ar:isobutane:CF4 (95:2:3), with a record diffusion constant as 
low as 26 //m/i/cm, the situation is even more catastrophic: the hodoscope effect is felt over all 
distances even with 1 mm pads (Fig. 12). 

This study shows that, to obtain the target resolution required to fulfill the ILC physics program, 
either smaller, digital pads are necessary, or a spreading of the charge onto several pads has to be 
implemented after amplification. The latter can be achieved either in a multi-GEM structure by 
maximizing the natural defocusing in the transfer between two GEMs, or by the use of a resistive- 
capacitive continuous layer (resistive foil) as proposed by some of us [14] . 



5 Conclusions 

Successful operation of a Micromegas TPC in a test beam enabled us to measure its spatial resolution 
under a magnetic field as well as the properties of an Ar isobutane mixture. The drift velocity and 
the transverse diffusion constant for at S = 0, 0.5 and 1 T are found to be in good agreement with 
the predictions of Magboltz. The obtained spatial resolution as a function of drift distance is well 
reproduced by an analytical calculation developed in this paper. The effective number of electrons 
(A'eff) is measured to be 18.5 ± 1.1, which is consistent with an estimate based on the primary 
ionisation statistics and the gas gain fiuctuations. 

The extrapolation to the conditions of the Linear GoUider shows that the goal resolution of 100 
microns cannot be reached with 2.3 mm pads with the technique used in this test, and that another 
technique will probably be necessary (digital smaller pad readout or resistive anode readout, for 
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Figure 12: Resolution in the Linear Collider case as a function of the drift length, for two values of 
the pitch, for a fast gas mixture with a large collision time. 

instance) . 
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A Pad Response Function in the Large cTd limit 

From Eqs. (15) and (17) we obtain the average charge on pad j: 

{Qj{x)) = NG {fj{x + Ax)), 
resulting in the average charge fraction on pad j: 

{Qj{x))/{NG) = {fj(x + Ax)) = ldAxPD{Ax;ad)fj{x + Ax) 

( 1 /A,T^ 2' 



fwj—x+w/2 ][ 

dAx—i= — exp 

lwj-x-w/2 \/lT^(7d \ 2 \ arf 



rwj—x- 
Jwj—x — 

r+w/2 

/ . 

J-w/2 



^ex ( + 



In our standard pad response function analysis, we plot this as a function of the pad center as 
measured from the average charge centroid: 

x{x) = jw - {x{x)) = jw-Y, (kw) {fk{x + Ax)) 

k 
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jw — X + O 



w 



~ JW — X. 



In the large limit, the charge fraction hence has the following functional form: 

2\ 



1 /-Wa 1 

W J-w/2 V271 



exp 



1 /,x + e 



Since Qpr{x) is apparently normalized to unity, we have 



(30) 



-)x^ 



(^PR = / dxQpR{x) 

I r+w/2 r+oo \ 



W J-w 2 J-i 



dx r— — exp ( — - f — ) I x^ 



-w/2 J-oo \/27r(Td ^ \ 2 \ ad 

W J-w/2 12 



and, consequently. 



12 



(31) 



B Scaling and (t^ ^ oo Limit 



As long as the avalanche can be regarded as point-like, and the pad response function is hodoscope- 
hke as given by Eq.(23), any observable with the dimension of length should scale as aa times a 
function of {w/ad) or w times a function of (w/aa). In the case of the first term of Eq.(25), it should 
hence scale as 



-1/2 \W 



(32) 



where the pad pitch (w) can only appear in the function through the ratio: crd/w. This term 
represents the well known 5'-shape systematic bias in the charge centroid for a finite pad pitch and 
will vanish in the w — > limit: Fi(oo) = 0. It is, however, non-trivial whether this will vanish in the 
(7(i — > oo limit for a fixed pad pitch: w — constant. We can show analytically that this is indeed the 
follows: 



J(j-l/2)w-x J{-j-l/2)w-x 
■ i'{j+l/2)w-x p{j+l/2)w+x 

J(i-l/2)w-x J(: 



(-j+l/2)w-x' 



r(j+l/2)w-x 

Y,ijw){fj{x + Ax)) = / dAxPn{Ax;a,) 

j j J(j-l/2)w-x 

" l-(j+l/2)w-x 

oc 

oo 
3=1 



roo roo 
J{j-l/2)w-x J{j-l/2)w+x 



{j-l/2)w+x 



dAxPD{Ax;ad) 
dAxPoiAx-Gd) 



dAxPD{Ax;ad) 
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r{j-l/2)w+x 



7=1 



l/2)w-x 



dAx PD{Ax;ad), 



where we have used the fact that Pd{Ax; aa) is an even function of Ax. 

Noting that —w/2 < x < -\-w/2 and hence \x\ <^ a a in the large a a hmit, we can further the 



calculation by Taylor expansion: 



r\3- 

Y^ijw) {fj{x + Ax)) - / dAxPD{Ax;ad) 



~ 2x'£wP„((j-l/2)w;ai) 

~ / rf^PD((j-i/2V + e;f7d) 

/•oo 

= 2x / ciAxP£,(Ax;(Jd) = x. 
Jo 



In order to see the large aa limit of the second and the third terms of Eq.(25), we need to evaluate 
the following: 



I = E(j^)' + Ax))- 

n(j+l/2) W-X 



x' 



w 



0-1/2) 



dAxPoiAx; aa) - x^ 



- E 



w/2 
+w/2 



dC (jw + ifPoiiw + i-x;aa)- 

-E / di 2{jw)i + e PD{jw + ^-x;ad) 
= J dxx^ Pd{x — x; ad) — 



- E Pd{3^ - x; aa) / d^ 2{jw) ^ + ^ 



+w/2 



w/2 



- - E PdUw - x; Gd) [2 ({3w)x - {jwf) + cr; 
j 

...2 

E [2 {{3w)x - {jwf) + cr^] Poijw - x; aa) w 



— a 
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r / NT r+w/2 



w 



Viol 



[2x^ -21 + al 



w 



21 -at. 



where use has been made of Taylor expansion to the first order of ^ : 



(33) 



Poijw + ^-x;ad) ~ PdUw - x; aa) - Poijw - x; a^) 
= Pd^w - x\ad)\l- ''^ 2^ i] 



JW — X 



(Ja 



and the fact that the odd functions of ^ vanishes upon integration and x^ < Substituting the 

right hand side of Eq.(33) in Eq.(33) iteratively and ignoring the terms of 0((w^/(tJ)^), we finally 
arrive at 



7 ~ aj + 



w 



12^1 



ICTa — (Ta 



w 
12 



(34) 



Notice that cirf-independent term can be regarded as proportional to {a^ x (w/cr^))^, and hence 
is no exception for the scaling law. 



C Geometric Mean Method 

In this appendix we give a simple demonstration of the geometric mean method applied in the 
analysis to estimate the spatial rcsohition. 

First, in the case where the hit point in question (say, the i-th point Xi) is excluded in the track 
fitting, the residual is given by 

^Xi Xi Xi , 

where Xi represents the estimator for the i-th point given by the track fitting using the remaining 
hit points. Its variance is 

aL, = ((Ax,)') = 4 + 4, (35) 

the sum of the true spatial resolution and the tracking error. 

Next, in the case where the hit point in question is included in the track fitting, the estimator for 
the i-th hit point is given by the weighted mean of Xi and xf. 

Xi , 

with w^. {wxj being the corresponding weight: l/cr|. (1/cr^.). The residual is hence given by 

2 2 

A / ^/ ^Xi / ^ \ ^Xj A 

ZAXj = Xi Xi - ■ yXi Xij - ■ • L^Xi . 

Ot. + atL. + (7^2 

The variance of the residual in this case is therefore 

= {{Ax'i)') = {(xi - x'if) = ^^f^ . (36) 



3^ j 
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Combining Eq. (35) and Eq. (36), one gets the following relation 
which immediately gives the expression for the true spatial resolution: 

Cxi — V'^excl • CTincl ■ (37) 

It should be noted here that Eq. (37) is vahd not only for straight hne fitting as dealt with in the 
appendix of Ref. [15] but also for track fitting with any function such as circles used in our analysis. 
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